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ABSTRACT 

We present model fits to spectral energy distributions in the optical and NIR of > 100 
flat-spectrum radio quasars from the Parkes Half-Jansky Flat-spectrum Sample. We 
find that ~ 40% of the sources have power-law SEDs, while a similar number show 
evidence for two primary components: a blue power law and optical synchrotron emis- 
sion. The blue power law is similar to the dominant component observed in the spectra 
of optically-selected QSOs. There is strong evidence that the synchrotron component 
has a turnover in the UV-optical rest frame of the spectrum. In the remaining sources 
it is likely that the synchrotron peaks at longer wavelengths. This mixture of two 
components is supported by optical polarisation measurements in a subgroup of the 
sources. The sources with power law SEDs show evidence for an excess number of red 
power law slopes compared to optically-selected quasars. There are additional spectral 
components in some of the sources, such as dust and the underlying galaxy, which have 
not been considered here. 
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1 INTRODUCTION 

The spectral energy distributions of Active Galactic Nuclei 
provide a sum of the emission processes contributing to the 
energy output of the AGN. Clearly, several mechanisms con- 
tribute at most frequencies. If we can determine the different 
components as a function of frequency, then we can deter- 
mine the energy generation mechanisms which are impor- 
tant. 

Flat-spectrum radio quasars are known to be dominated 
in the radio spectrum by synchrotron emission. It is less clear 
how far these synchrotron spectra extend towards higher 
frequencies, or the magnitude of the contribution they make 
to the optical and near-infrared (NIR) spectrum - a region 
dominated by the Big Blue Bump. 

The far-IR emission of flat-spectrum quasars appears t o 
be dominated by synchrotron emission (Haas et al. 1998), 
consistent with a Doppler-boosted component that swamps 
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NIR wavelengths comes from both the variability in the NIR, 
and from the strong polarisation of a few sources. Stickel et 
al. ( 199(^ ) found that variability was a general characteristic 
of optically faint flat-spectrum radio sources, although they 
did not specifically discuss synchrotron emission. Instead, 
they suggest that reddening due to intervening galaxies or 
the host galaxy itself was the likely cause of the red optical 
colours. 

While radio jets are quite common, and have been 
widely observed, optical counterparts to these jets are quite 
rare. In fact, only 14 optical jets are currently known (see 
O'Dea et al. (1999) and references therein). The features (i.e. 
emission knots) in the optical jets match the positions of fea- 
tures in the radio jets, indicating that the emission regions 
and mechanisms for the two spectral regimes are associated. 
Optical jets can be used to put strong constraints on the 
particle energetics, by providing limits on the maximum en- 
ergies, and on the acceleration mechanisms, due to the short 



emissiofi in the optical ha.s been considered for many years 
to be necessary to explain the optical to NIR photometry of 
optically faint and red flat-spectrum radio quasars. This has 
mostly b een on the basis of the steepness of t he optical-NIR 



spectra (Rckc et al. 1979; Beichmann 1981 Bregman 1981 



Rieke et al. 1982). The steepness of these spectra suggested a 
cutoff, or sharp break, in the electron energy spectrum. Sup- 
porting evidence for the presence of synchrotron emission at 



lifetim es of particles at these energies ( Meisenheimer et al 
1996) ). 

The quasars considered in this paper are selected 
f rom the Parkes Half-Jan sky Flat-spectrum Sample (PHFS) 
(Drinkwater et al. 1997), which consists of 323 objects se- 
lected to be radio-loud (S2.7 GHz > 0.5 Jy), and have a flat 
radio spectrum (a ? 7/5.0 > —0.5, S(u) o c v a )- These quasar s 



have been shown ( Webster et al. 1995] ; Francis et al. 200C) 
to have a large spread in B — K colours, with the reddest 
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objects having B — K> 7. Masci, Webster & Francis (199£ 



showed that this spread could not be accounted for by the 
emission from the host-galaxy. 

In this paper we test the idea that optical synchrotron 
emission causes the red optical-NIR colours of the PHFS 
quasars. Models representing emission from both an opti- 
cal synchrotron component and a blue optical power law 
(representing the continuum emission from an unreddened 
quasar) are fitted to broad-band optical-NIR spectra. Stan- 
dard goodness-of-fit techniques are used to determine mod- 
els which are consistent with the observations. 

The d ataset was compiled by Francis, Whiting & Web- 
ster (200C) (FWW hereafter), and comprises broadband op- 
tical and NIR photometry in the bands B, V, R, I, J, H and 
K. This photometry is quasi-simultaneous, meaning all ob- 
servations for a given source were made within several days 
(six at most) of one another. This minimises the effects of 
source variability. The source selection used for this paper 
is explained in Section ^. The simultaneity, as well as the 
breadth and density of the spectral coverage of this dataset 
provides an excellent basis to model the broadband emission 
from a large number of flat-spectrum quasars. 

The dataset is described in Section ^| In Section ^ we 
test the hypothesis that the emission in the optical-NIR re- 
gion is well-fitted by a single power law. In Section ^| we 
describe the more sophisticated model that was fitted to the 
data - representing the accretion disk and synchrotron emis- 
sion - and the method that was used to perform the fitting. 
The results of this fitting are described in Sections ^ & j| 
while further tests of the model fits are examined in Sec- 
tion ^, using polarisation and emission line measurements. 
A possible rival model - that of black-body emission by hot 
dust - is considered in Section ^. The effect of emission lines 
on the photometry and the resulting fits is investigated in 
Section ^. Finally, Section [l^ contains discussion of the re- 
sults and their implications. 



2 DATA 

The data which have been fitted by these models are de- 
scribed in detail in FWW. A total of 157 sources from the 
PHFS were observed, with quasi-simultaneous broadband 
photometry observations in the bands B, V, R, I, J, H and 
K. These magnitudes were converted into broadband fluxes 
using the zero points given in the same paper. 

2.1 Excluded sources 



Not all sources in FWW were used. Firstly, we only con- 
sidered those sources with complete contemporaneous pho- 
tometry in all bands, or those with only one observation 
missing (either not observed or an upper limit). This was 
to ensure that the number of degrees of freedom in the 
model fitting was greater than zero. Those sources with more 
than one band missing were not included in the analysis. 
Sources without a measured redshift were also excluded, 
since the redshift is needed to obtain the correct shape 
of the observed synchrotron spectrum. One of the sources 
t hat h ad an unknown redshift (z — 0) in Drinkwater et al. 
( h997|), 0829+ 046, has a published redshift value of z = 0.18 
(Falomo 1991), which is used here. Also excluded were low 



redshift galaxies, which had a prominent 4000A break be- 
tween the B and V bands. These sources show strong ev- 
i dence for flux from the underlying galaxy in their spectra 
(Masci et al. 1998), and thus need an additional galactic 
component to be modelled accurately. Sources with z > 3 
were also excluded, as these had strong Lya breaks present 
between B and V. This reduced the number of sources from 
157 to 117. 

Twenty-one of the reddest sources (B — K > 5) in the 
subsample were amongst those excluded, as they have only 
upper limits on B, V and possibly R, or no redshift. The 
sources with upper limits are both red and optically faint. 
Since the flux of these sources typically decreases rapidly in 
the blue (i.e. from V to B), there are three possible expla- 
nations for them. They are either dust-reddened (causing 
the blue decrease in flux), high redshift (and hence Lya ab- 
sorbed), or they are dominated by synchrotron that turns 
over rapidly. Future analysis of the spectral energy distribu- 
tions (SEDs) of these sources will only strengthen our final 
conclusions. 

Those sources excluded purely because of their lack of 
measured redshift generally have high B — K colours (B — 
K > 5), since the most likely reason they have no redshift 
is that they are faint in the optical (and in particular in the 
blue). Further discussion of the sources excluded from the 
sample is given in Section 



2.2 Errors on photometry 

The photometry given in FWW quoted error bars, where the 
estimated error comprised two parts: a random error com- 
ponent and an assumed 5% error in the photometric zero 
points, which were added together in quadrature. The pho- 
tometric zero point errors were estimated from the scatter 
in zero points between different standard star measurements 
in an individual night: Francis et al. adopted a value of 5% 
to account for this error. 

However, this zero point error ignores a number of fac- 
tors that we believe may be important for the analysis in 
this paper. These factors are as follows: 

• The optical zero magni tude fluxes were taken from 
Bessell, Castelli & Plez (1998). These fluxes were derived for 
an AO star, and so will be slightly incorrect for a quasar spec- 
trum. This has the effect of introducing small colour terms 
into the photometry, the size of which will depend on the 
spectral index of the object being observed. A similar effect 
will of co urse be present in the NIR. Bersanelli, Bouchet & 
Falomo (1991) found that spectral shape differences could 



produce systematic errors of at least a few percent in the 
NIR flux. 

• The zero point f luxes in the optica l are taken from 
a different reference (Bessell et al. 1995) to those in the 
NIR, which were calculated by P.McGregor (CASPIR man- 
ual, MSSSO, ANU), assuming that Vega is well represented 
in the NIR by a black body of temperature 11200 K, and 
normalisation F x f 555 nm) = 3.44 x 10~ 12 W cm -2 /im _1 
( Bersanelli et al. 1991 ). These different zero points may not 
be exactly equivalent, which will produce small offsets be- 
tween the optical and NIR parts of the SED. 

• If the sky conditions at Siding Spring Observatory were 
not completely photometric for all observations for a given 
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source (particularly if the transparency changed between dif- 
ferent bands) , then the measured photometry will have small 
band-to-band errors present. 

• The observations for each source were taken quasi- 
simultaneously (meaning all observations were made within 
at most a six-day period) , to minimise the effects of variabil- 
ity. However, a number of the sources in the PHFS h ave been 
found t o exhibit intra-day variabil ity in the optical ( Romerc 
et al. l!)99; Hcidt fc Wagner 1996| ), and a larger number no 
doubt have similar properties to these. Therefore, variability 
on timescales of the order of those separating our observa- 
tions is likely for some of the sources. Such variability can 
be up to 0.1 mag over the period of a night. 

• Finally, the presence of strong emission lines in the 
quasars' spectra could boost the flux of a band above the 
level of the continuum. We have tried to quantify this effect 
in Section [], although this is not possible to do for many 
sources, due to the lack of both photometry and a spectrum. 
Some additional effects may be due to line blends (such as 
Fell blends), but this is probably not as strong an effect as 
you would find in a radio-quiet sample. 

Hence, to take account of all these factors, we have in- 
creased the systematic error in the photometry to 10%. The 
random error is kept at the same level as that presented in 
FWW, and is added in quadrature to the systematic error. 



3 POWER LAW MODEL 

The object of our analysis is to find physical models to ex- 
plain the optical and NIR emission. FWW found that about 
90% of the PHFS have approximately power law SEDs. We 
first wish to test this more rigorously. As a starting point, we 
choose to fit (naively) a simple power law, with an uncon- 
straine d spectral index. This will separate out the sources 



that have power law SEDs from those that show some cur- 



vature in their spectrum. 



3.1 Fits to data 

The model we choose to fit is /pl(A) = c\ a , so that the 
normalisation c and the spectral index a vary. (Note that 
this model implies f v oc v~ 2 ~ a .) 

This model is fit to the data using a least-squares 
method. This generates a \ 2 value: 



X 



E 



(Vi - fpL(Xi)) 2 



which indicates the goodness of fit. A fit to a source will be 
deemed to be "good" when the value of \ 2 m l ess than the 
cut-off value corresponding to the 99% confidence level of 
the x 2 distribution. For a source with 5 degrees of freedom 
(as is the case for most sources with this power law model), 
this cutoff level is 15.09. (Note that increasing the value of 
this cutoff is equivalent to increasing the confidence level 
- for example, the 99.5% cutoff is 16.75 for 5 degrees of 
freedom.) If the \ 2 value is greater than this cutoff level, 
then we reject the null hypothesis that the power law model 
fits the data. 

When we fit this power law model to the data, we find 
that 83 sources (or 71% of the total) have good fits. The 




Opticol index 

Figure 1. Histogram of fitted power law indices. Hatched his- 
togram indicates good fits, while open histogram indicates all 
fits. 



distribution of resulting power law indices, both for the good 
fits and for all sources, is shown in Fig. [I]. 

The spectral indices of these good fits span a wide range 
of values. At one extreme there are the sources with rela- 
tively blue SEDs (a < —1.4). These sources are charac- 
terised by their blue continuum, the presence of moderate 
to strong emission lines, and ge nerally low X-ray f lux (most 
were not detected by ROSAT ( [Siebert et al. 1998| )). 

At the other extreme are the sources with redder SEDs 
(that is, flatter in fx), with — 1 < a < 0. These sources are 



blazar ty pe objects, with high optical polarisation (Wills et 



al. 1992) as well as relatively weak (or even absent) emission 



lines - in fact, all the BL Lac objects in our sample are in 
this region. 



3.2 Interpretation of power laws 

So, we have fit a power law to a large majority of the sources 
in our sample, spanning a wide range of spectral indices. 
Are the physical processes that generate this power law the 
same for all sources? That is, does the power law in the blue 
sources have the same origin as that in the red sources? 

The first class of sources - the blue sources - are being 
fit by a blue power law, which has similar colours to the 
blue power la w emission seen in optically-selected quasars 
(Francis 1996). This is likely to be the optical part of the 
continuum emission from the accretion disk (often termed 
the big blue bump). 

However, the power law being fit to the redder sources is 
most likely of different origin to that seen in the blue sources. 
These objects exhibit characteristics commonly associated 
with optical synchrotron emission (such as high optical po- 
larisation and lack of prominent emission lines), and so we 
postulate that this emission is, at least in part, some form of 
synchrotron emission. The slope of the power law can then 
be used to determine p, the power law index of the electron 
energy distribution (i.e. defined such that N(E) oc E~ p : see 
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Section |4.2| ). Using the values shown in Fig. |l], we obtain 
2 < p < 6 (using —1.5 < a < 0.5 and p = 2a + 5). 

If synchrotron emission is present in the spectra of at 
least some of these quasars, then we can ask the question 
"Is the synchrotron component best modelled by this power 
law?" The synchrotron component will be present in one of 
two forms: a power law caused by an unbroken (power law) 
electron energy distribution, or a turning-over component 
caused by a break or a cut-off in the electron energy distri- 
bution. (Note that the power law can also be produced by 
a synchrotron spectrum turning over at higher frequencies 
than those observed.) 

Both of these forms can be tested. The power law model 
must produce power law indices that are consistent with 
slopes of plausible energy distributions. The presence or 
otherwise of a turnover can be evaluated by examining the 
sharpest possible turnover (caused by an abrupt cut-off in 
the energy distribution at some maximum energy). This will 
provide the maximum contrast with the power la w, and is 



consist ent with m odeling done by other authors (Meisen- 
heimer et al. 1996, for example). 



Many of the sources, while they have power law fits 
that we can not reject at the 99% confidence level, show 
evidence for curvature in their SEDs. This curvature can be 
'n'-shaped (higher in the middle than at either end), 're- 
shaped (lower in the middle) or perhaps take the form of 
an inflection (e.g. the flux decreases, levels off and decreases 
again). (See Fig. ^ for illustrations of the different types.) 
Most of the sources classed as BL Lac objects are 'n'-shaped, 
and so we postulate that these sources are dominated by a 
synchrotron component that is turning over in the optical. 

Other sources, however, are bluer in the optical than in 
the NIR (i.e. show an inflection, without the turn-up seen 
in 'u'-shaped sources). This is a possible indication of the 
presence of excess emission in the NIR, in addition to a blue 
power law. We propose that this excess emission is due to 
a synchrotron component that has turned over in the NIR, 
and so does not dominate in the optical, where the dominant 
emission is instead a blue power law similar to that observed 
in the bluest sources. 



Over the wavelength of our observations, the BBB is 
modelled as a simple power law, f\ oc X aB or /„ oc u~ 2 ~ aB . 
Francis (1996) found the median slope for a subsample of 
quasars from the Large Bright Quasar Survey (LBQS), taken 
from optical/NIR photometry, was /„ oc j,- ' 35 * - 3 ^ anc [ 
noted that the observations were consistent with an intrinsic 
continuum slope of /„ oc v~ ' 3 that is reddened by various 
amounts of dust. We therefore take our value of olb to be 
— 1.7. The effects of allowing the value of as to vary are 
considered in Section [l^. 



4.2 Synchrotron emission 

All the PHFS sources are radio-loud flat-spectrum sources, 
and thus very likely have relativistic jets that emit syn- 
chrotron radiation, at least at radio frequencies. Could this 
synchrotron emission extend up to the optical/NIR part of 
the spectrum? Our power law fitting from the previous sec- 
tion provides circumstantial evidence for this: the redder 
sources in the optical tend to be the ones with higher po- 
larisation (a good sign of synchrotron emission) and less 
prominent emission lines (possibly a sign that the emission 
lines are being swamped by the presence of a synchrotron 
component). These pieces of evidence are investigated more 
deeply in Section ^| 

As discussed above, the synchrotron spectrum could 
take the form of either a power law, from a power law dis- 
tribution of electron energies, or a power law with a break 
or turn-over, due to an electron energy distribution that ex- 
hibits a break or even a cut-off. This latte r type of spectrum 
has been seen in o ptical synchrotron jets ( O'Dea et al. 199£ ; 



Scarpa et al. 1999] ) , where the optical spectrum is like 



7H3 



compared to a radio-optical spectrum of v 



Additionally, a synchrotron spectrum that has a turn- 
over will, when combined with the blue power law, be able 
to reproduce an inflection-like SED. Such a spectrum, par- 
ticularly in the region of the turn-over, will also be quite red, 
thus accounting for the red colours of many of the SEDs. 



4 PHYSICAL MODELS 

In response to this phenomenological classification, we pro- 
pose the following physical model. There are two compo- 
nents in this model: one is a blue power law, representing 
continuum emission from the accretion disk region; and the 
second is synchrotron emission, representing emission from 
the relativistic jet (that we know to be present due to the 
flat-spectrum radio emission seen in these objects). 



4.1 Accretion Disk emission 

We find, from the simple power law fitting, that the bluest 
sources have power law continua. The slopes of these power 
laws are consistent with them being the same component as 
that seen in optically-selected quasars, commonly termed 
the "Big Blue Bump" (BBB). We take this component 
to represent the underlying quasar continuum in the UV- 
optical part of the rest-frame spectrum - that is, the "un- 
reddened" quasar continuum. 



4-2.1 Analytic modelling 

We consider here synchrotron emission from a population 
of electrons with an energy distribution with the form of a 
power law up to some maximum energy and zero beyond 
this (i.e. an energy spectrum with an abrupt cutoff). This 
can be expressed as a distribution of the Lorentz factor 7 of 
the radiating electrons: 



N("f) d-y oc 



K7 v d7 1 < 7 < 7 C 







7 > lc 



T his is consistent with modelling done by Meisenheimer 
et al. (Il99d) on the jet of M87. They found that the overall 



synchrotron spectrum of the brightest parts of the jet was 
best described by a spectrum that had a sharp cutoff at 
v c w 10 15 Hz, with an energy distribution of the form of a 
straight power law N(~y) oc 7~ 2,31 , with a rather abrupt high 
energy cutoff. We consider the effect of using a power law 
synchrotron spectrum instead in Section [l^. 

Such a synchrotron spectrum is straightforward to 
model analytically. We use the "classical" synchrotron 
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model, as first calculated by Schwinger (194£), and as de - 
rived by a nu mber o f authors, particularly Pacholczyk ( 197C ) 
and Longair (1994), whose derivation we follow. 

The single particle luminosity for a radiating electron 
(mass m e , charge e, Lorentz factor 7 = E/m e c 2 , and pitch 
angle, or angle between trajectory and magnetic field direc- 
tion, 9) in a uniform magnetic field B is given by 



Lip(w) = 



^3e 3 Bsin0 



F(x) 



8iv 2 eQcm e 

where u) = is the angular frequency, 

2uj/3m e 
37 2 e_B sin 8 ' 

and F{x) is denned in terms of the Bessel function ^5/3(2) 
by 



F(x)=x K 5/3 (z)dz. 

J x 

We are interested in the luminosity of a population of 
particles, so we need to integrate L\ p {ui) over suitable distri- 
butions of energies and pitch angles. The energy distribution 
is that given above, while the pitch angle distribution that 
we use is an isotropic one, where the probability distribution 
is p(6)dd = i sin6> d8. Thus, the integrated luminosity from 
such a population is 



167r 2 eocm e 



7" p F(s) d7 d8 



An example of such a spectrum is shown in Fig. ^, for 
7 C = 10 4 and B = 1(T 4 T = 1G (the value of k has been 
taken to be 1). The peak frequency v c depends on these two 



valu es, and can be shown by simple arguments (Blandford 
199C |) t[ a be approximated by v c ~ "/ 2 B MHz (where B is 
measured in Gauss). 

The slope of the power law tail (at frequencies v <C v c ) 
is related to the energy power law index by as = (p — 5)/2 
(where f oc X as ). The energy distribution for the spectrum 
in Fig. H is taken to be N(j) oc 7~ 2,5 (i.e. p = 2.5), giving a 
power law of f\ oc A -1,25 . 



We consider here a range of p values from p — 2.0 to 
p = 3.0, which gives a range of long-wavelength power law 
slopes of as = —0.5 to as = —1.0. This range covers the dis- 
tribution of r adio-to-optical slopes observed in optical syn- 
chrotron jets (Scarpa & Urry 1999). Allowing p to vary does 
not significantly alter the results of our analysis - see Sec- 
tion hd for further discussion. 

We also note here that a value of p > 3 means that 
the vF v flux will increase towards longer wavelengths (since 
F v oc z,-^" 1 )/ 2 and so vF u oc z,"^- 3 )/ 2 ) and this results 
in the radio flux being severely overestimated by the fitted 
synchrotron component, since the radio emission always has 
a lower vFv flux than the optical. (This assumes that the 
same synchrotron component is responsible for both the op- 
tical and radio emission, which is an assumption commonly 
mad e, particularly for the mod elling of optical synchrotron 
jets ( Meisenheimer et al. 1996).) 



4.3 Model fitting 

These two components (the blue power law and the syn- 
chrotron component) are combined linearly to form a model 
/c(A) = a\~ 17 + bf Bync h{X) that is fit to the data in the 
same way as the power law model (that is, using \ 2 minimi- 
sation). The reduced x 2 value (that is, X 2 / v ) f° r each of the 
two models (combined and power law) are compared, and 
the model with the lowest x 2 l v i s chosen to be the best fit 
model. 

In fitting the combined model, the location of the peak 
wavelength of the synchrotron spectrum, A p , was allowed to 
vary. This variation was allowed to occur over a range of 
rest frame wavelengths such that the curvature of the spec- 
trum caused by the turn-over affected the synchrotron flux 
in the region of the data points (in other words, we did not 
want to just be fitting the power law part of the synchrotron 
spectrum) . Quantitatively, we took the minimum peak wave- 
length to be half a decade shorter than the B band (0.44/im) 
shifted to the rest frame. We then considered 20 A p values 
per decade (evenly spaced in log 10 A p ), up to a maximum 
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ALL SOURCES 



POWER LAW 



n - COMBINED 




Power low index a 

Figure 5. Distributions of the fitted power index a. The open 
histogram shows the results of the power law fit to all 117 sources, 
the hatched histogram shows which sources have a power law fit 
accepted at 99% confidence, while the cross-hatched histogram 
shows which of those sources have power law fits better than the 
combined fit. 



Figure 4. Histograms of the reduced \ 2 values, for all sources 
and for each of the best-fit models. The dotted histogram shows 
those sources whose fits are rejected at the 99% confidence level. 



peak of 10/xm. For each of these synchrotron functions, a 
best fit to the data was found, and then the best of these 
was chosen, giving the best fit \ p value for that source. 



While many sources have been fitted better with the 
combined model, a large number are still preferentially fit 
with the power law. If we plot a histogram (Fig. ^) of the 
power law indices of those still fit by the power law model, we 
can see that those sources that are preferentially fit by the 
power law model are the bluer sources, while the majority 
of the sources with indices a > — 1 are fit better by the 
combined model. 



5 RESULTS FOR PHYSICAL MODELS 



6 FITTED SYNCHROTRON COMPONENTS 



This combined model was fitted to the photometry, and com- 
pared to the power law fits. The best fitting model was cho- 
sen on the basis of the lowest reduced \ 2 value, as described 
above. For the default values of the parameters (p = 2.5 and 
ub = —1-7), we find that 93 sources (or 79% of the total) 
are well fit by one of the models. Of these, 48 are best fit 
by the power law model, and 45 by the combined model. 
How these numbers c hange with different parameter values 



is discussed in Section 10.4. A selection of the fits are shown 
in Fig. ^, for a range of power law slopes and synchrotron 
peak wavelengths. 

We note here that although there are 48 sources best 
fit by the power law model, many of these have combined 
fits that are only slightly worse than the power law fit. This 
indicates that there is not a great deal of difference between 
the fits of the two models. This is not the case with the 
combined model sources, as for most of these the combined 
model fit is a lot better than the power law fit. 

A histogram of \ 2 l v values is shown in Fig. ^, with the 
distributions for the two different models shown separately. 
The distribution for the sources best fit by the combined 
model is noticeably broader than that for the power law 
sources, with more sources having very low x 2 values. 



The properties of the synchrotron components that are fitted 
as part of the combined model are of particular interest. As 
can be seen in Fig ^, the peak wavelengths are restricted to a 
relatively narrow range of values (approximately a decade in 
wavelength). However, this is likely to be largely a reflection 
of the distribution of the wavelengths of the photometric 
points. 

The strength of the fitted synchrotron component varies 
considerably from source to source. In Fig. |^, the ratio of 
the synchrotron and power law components at a rest-frame 
wavelength of 0.5/im is shown for all sources best fit by the 
combined model. The main bulk of this distribution spans 
nearly four orders of magnitude. This large range of values, 
which is also seen in the normalisations of the individual 
components, indicates that we are seeing a continuum of 
variations of these components, probably due to variations in 
the strengths of the inner jet and emission from the accretion 
disk and/or surrounding regions. 

We also note that a small number of the sources at the 
high-ratio end of the distribution are faint, red sources, that 
are likely to be significantly dust-reddened. They are thus fit 
with a dominant synchrotron component, as the synchrotron 
spectrum has the approximate form of a power law with an 
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Figure 3. Examples of the various types of photometry, and the fits to them. The horizontal axis on each plot is wavelength, in /im, and 
the vertical axis is fx in W m -2 fixa~ . The locations of notable emission lines are shown, based on the redshift of each quasar. For the 
combined fits, the individual components - the power law (ctg = —1.7) and the synchrotron component (p = 2.5) - are shown as dotted 
lines. The An v alue, where given, is for the quasar rest frame. A full compilation of the fits to all the quasars in the sample is given in 
Whiting (hood). 



exponential cut-off (which is the same as a power law with 
dust extinction). 



7 TESTING THE FITS 
7.1 Polarisation 

We have shown that a large fraction of the sources in the 
PHFS show evidence for the presence of optical synchrotron 
emission, where the amount of synchrotron emission present 
in the spectrum changes with wavelength. How else can we 
test this model? One of the key features of synchrotron radi- 



ation is its high degree of polarisation. If there is significant 
synchrotron emission at optical and NIR wavelengths, then 
one would expect to be able to detect a corresponding po- 
larisation. Indeed, this has been used as a way to confirm 
t he pr esence of synchrotron emission in optical jets (Baade 
( |l956| ) provided the first example of this for the jet of M87.) 



In our combined model, the synchrotron component is 
the only polarised component, as we assume that the BBB 
component, which is essentially emission from the accreti on 
disk, is unpolarised (P < 1% for the BBB ( [Antonucci 198^ )). 
Thus, the amount of polarisation will depend on the propor- 
tion of the total flux that is due to the synchrotron emission. 
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A P eak (Observed) 

Figure 6. Histogram of peak wavelengths, showing all sources fit 
best with the combined model, in a) the rest frame, and b) the 
observed frame. The solid histograms show all sources accepted 
at 99% confidence, while the hatched histograms show the BL 
Lac objects. The locations of three prominent emission lines and 
photometric bands are given for reference. 




10 1 1 10' 
0.5 jOtm Ratio 

Figure 7. Histogram of ratios of the synchrotron component 
to the power law component of the combined model, for those 
sources best fit by the combined model. The ratio is calculated at 
a rest-frame wavelength of 0.5^tm. The hatched histogram shows 
the BL Lac objects. 



• Combined model 

* Power law, a>— 1 .3 



0.5 1 
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Figure 8. Polarisation from Wills et al. (1992) as a function of 



the proportion of the total flux made up by synchrotron. The 
sources are given different symbols according to the nature of 
their best fit model. 



Furthermore, if the relative amount of synchrotron emission 
changes with wavelength, as it does with the models we have 
fitted, then the amount of polarisation should also change 
with wavelength. 

Such dependencies have been investigated previously by 
a number of different authors, for sam ples th at include some 
sources considered here. Wills et al. (1992) studied a large 



sample of bright, flat-spectrum core-dominant quasars, mea- 
suring their optical polarisation. An interesting result is that 
they found that the fraction of quasars with P > 3% in a 
fixed observed passband decreased with increasing z, pos- 
sibly indicating that the percentage polarisation decreases 
towards shorter rest frame wavelengths. This would be con- 
sistent with the presence of a synchrotron component turn- 
ing over in the optic al reg ion. 

Impey & Tapia ( 1990 ) present radio and optical data for 
a slightly larger sample of radio-selected quasars, including 
optical polarisation measurements. They find strong statis- 
tical links between strong optical polarisation and properties 
such as compact radio structure, superluminal motion and 
weak emission lines. They explain this by requiring the op- 
tical emission, as well as the compact radio emission, to be 
relativistically be amed . 

Smith et al. ( Jl988| ) obtained multicolour (UBVRI) po- 
larisation measurements of 11 highly polarised quasars, and 
found that three of these exhibited decreasing polarisation 
toward shorter wavelengths, which they modelled as a com- 
bination of polarised synchrotron emission, and two unpo- 
larised components, from the broad-line region and the ac- 
cretion disk. None of these sources, however, are part of our 
sample. 



7.1.1 Optical polarisation 

Firstly, we wish to compare our model predictions with pub- 
lished optical polarisation measurements. We use th e larg e 
catalogue of measurements compiled by Wills et al. (1992). 
We want to compare these measurements with the fractional 
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amount of flux due to the fitted synchrotron component. 
However, since the observations of Wills et al. were made 
without a filter, we have determined the average synchrotron 
fraction by integrating over the range 0.3/^m - lfim. 

We have then plotted in Fig. | the percentage polarisa- 
tion as a function of this average synchrotron fraction. For 
those sources that were best fit by the power law model, 
we have calculated the fraction from the combined model 
fit, and indicated these sources by a different symbol. All 
the power law sources with polarisation measurements had 
relatively red power law indices (i.e. a > —1.3). 

The spread in polarisation measurements at large syn- 
chrotron fractions is much greater than at low fractions, in- 
dicating that the high polarisation sources generally have 
large amounts of synchrotron fitted to them (at the wave- 
lengths at which the polarisation is measured) . Additionally, 
all but one of the sources at zero (or near-zero) synchrotron 
fraction have low optical polarisation. 

The two exceptions to this picture are 0202—172 and 
1020-103. Firstly, 1020-103 has a very high synchrotron 
fraction (~ 95%), but has very little optical polarisation 
(P = 0.58%). In this case, the source has a power law con- 
tinuum with a slight curvature, which is fit well by an almost 
pure synchrotron curve. It may be that this curvature is due 
to other effects than synchrotron, which explains the lack of 
polarisation. A possible candidate is contamination from the 
very strong Haline, which would boost the continuum level 
in the centre of the SED. 

Secondly, 0202—172 has a measured polarisation of 5.15, 
but also has a very blue power law continuum (a — —1.85). 
This continuum slope is possibly too steep to attribute to 
synchrotron (it implies a value of p = 1.3 - in turn implying 
a rather flat energy distribution ). However, the location of 
the synchrotron peak may have shifted in the period since 
the polarisation measurements were made (the polarisation 
measurements were taken 8 years prior to our photometry 
observations) . 

7.1.2 Near-infrared polarisation 

To try to avoid the problem of non-simultaneity of the pho- 
tometry and polarisation observations, we obtained polari- 
sation measurements of 8 quasars in the NIR, using the IRIS 
instrument on the Anglo Australian Telescope. The details 
and results of the observations will be presented elsewhere 
(Whiting et al., in preparation). These measurements are 
nearly simultaneous with the photometry measurements (a 
difference of ~ 40 days) , and so can be directly related to the 
fitted components. The wavelength dependence of the polar- 
isation can then provide an important test on the models we 
are fitting. 

Since the synchrotron component is the only polarised 
component, one would expect the percentage polarisation 
to be directly related to the amount of synchrotron flux 
present, and, in fact, the percentage polarisation will be di- 
rectly proportional to the ratio of synchrotron flux to total 
flux. In Fig. ^, we have plotted the polarisation of each of 
these quasars as a function of wavelength, and, on the same 
plot, the synchrotron ratio normalised (arbitrarily, as it is 
the wavelength dependence we are interested in, not the pre- 
cise normalisation) to the longest wavelength polarisation 
point (which is usually the 7f-band point). For complete- 



ness, we also show a line of constant polarisation, normalised 
to the AT-band point. All these sources are best fit with the 
combined model, with the exception of 1101—325, which is 
fit best by the power law model (its SED taking the form 
of a blue power law, A -1 ' 81 ). We have used the combined 
model fit to it for the purposes of Fig. |^. 

We have also included polarisation measurements from 
Wills et al. (1992) where they exist. Again, these are put 



at an observed wavelength of 0.5/im. Note, of course, that 
these points are not simultaneous with the NIR points. 

The fitted synchrotron component generally replicates 
well the wavelength dependence of the polarisation, al- 
though for some sources, such as 1313—333, the points are 
equally well given by a constant polarisation component. 
This is what you would expect from a pure synchrotron com- 
ponent, and these sources are typically BL Lac objects, from 
which you would expect to see a synchrotron-dominated 
SED. A notable exception to this is 0537—441, whose polar- 
isation is not fit well by a constant synchrotron component, 
but is fit better by a combination of a synchrotron and a 
significant power law component. However, it is apparent 
from these plots that having simultaneous optical polarisa- 
tion measurements would better help discriminate between 
the combined model and a constant polarisation model. 

7.2 Emission lines 

Synchrotron flux emitted from the PHFS quasars is most 
likely to come from a thin, relativistic jet, and the radiation 
will not be isotropic. This means that there will be very 
little synchrotron radiation directed towards the emission 
line clouds, which have a much larger covering angle. Thus, 
the ionisation of these clouds will be due to the continuum 
emission from the central accretion disk region - in other 
words, the Big Blue Bump. 

One might expect that adding a non-ionising syn- 
chrotron component to the continuum will have the effect 
of reducing the equivalent width of the emission lines from 
the BLR, due to the fact that the flux in the emission lines 
is not changed but the continuum flux is increased. 

To test this prediction, we compared the equivalent 
widths of five emission lines (Civ 1549, Cm] 1909, Mgn 
2798, H/3, and the doublet [O ill] 4959,5007) with the ratio 
of synchrotron to continuum flux at the line wavelength, to 
see if some form of an anti-correlation is presen t. The de- 
tails of t he observations will presented elsewhere (Francis et 
). Objects that had spectra taken were essentially 



al. 2001 



a random sample of the PHFS (subject to visibility during 
the observing run). 

In Fig. hoL we show the results of this comparison for the 
Mgn and H/3 lines, which are the two broad lines with the 
longest wavelength (and hence the two lines most likely to 
show a reduction in equivalent width). We have also plotted 
those sources best fit with the power law model. The value 
of the ratio used for these sources was taken from fitting the 
combined model, and so are upper limits to the ratio. 

The Mgn line does not show much relationship to the 
synchrotron ratio, while the H/3 line does show a reduction 
in equivalent width with increasing amount of synchrotron. 
This lends some support to the hypothesis that excess syn- 
chrotron light is present. The difference in the two plots is 
likely due to the presence of the turnover in the synchrotron 
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Figure 9. Near-infrared polarisation of a selection of quasars as a function of wavelength. Also shown are optical measurements from 
Wills et al. (1992). The solid lines show the fraction of the total flux made up by the synchrotron component (right-hand axis), while 
the dashed line shows this normalised (see text) to the fC-band data point, and is in polarisation units (left-hand axis). The dotted line 
shows constant polarisation, normalised to the i\-band data point as well. A * denotes the object is a BL Lac. 



flux, so that is has less effect at the shorter wavelength of 
Mgn. There are, however, several sources that have high 
synchrotron ratios together with high equivalent widths. The 
implications of these results are discussed in Section 



8 HOT DUST: AN ALTERNATIVE TO 
SYNCHROTRON? 

So far, we have shown that, for a number of sources, the 
optical-NIR photometry is well fit with a power law plus 
a curved component, which we have assumed to be the 
turnover of a synchrotron component. However could an- 



other model be used instead of synchrotron? One possible 
alternative is blackbody emission resulting from hot dust. 

To test this model, we used a blackbody emission spec- 
trum due to dust at a temperature of 1750 K (the sublima- 
tion temperature characteristic of dust grai ns co nsisting of 
graphite and silicates, e.g. Laor & Draine (1993)) emitted 



in the quasars' rest frame. A blackbody curve at this tem- 
perature would have its peak, in the quasar rest-frame, at 
1.66yjm. This blackbody spectrum was combined with the 
same A -1 ' 7 power law used in the combined model to pro- 
duce a model that was fit to the data. 

The fits generated by this model were almost always 
worse than those of the synchrotron model. This was due to 
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Figure 10. The equivalent width of a) the Mgn and b) H/3 lines as a function of the ratio of synchrotron to continuum flux at the 
emitted wavelengths. The sources are given different symbols according to the nature of their best fit model (the arrow indicates an 
upper limit for a source that was best fit with the combined model). Ratios for sources best fit by the power law model are calculated 
from the combined model fits, and so are upper limits to the ratios (for the crosses and stars). The dashed lines represent the expected 
change in equivalent width with increasing amount of synchrotron, for an emission line with intrinsic equivalent width of lOOA. 



the peak of the blackbody occurring at much longer observed 
wavelengths than the K band. In only 12 cases was the x 2 /u 
value from the dust model less than that of the synchrotron 
model, and in most cases this was because the dust model 
had one more degree of freedom. For each of these objects, 
the SED took the form of a blue power law (one object was 
A -1 ' 4 , and the rest were bluer than X~ ' ) that had a slight 
amount of reddening at the H and K bands, which was fit by 
the presence of the dust blackbody curve. For the majority 
of the sources, however, the hot dust model was a lot worse 
than the synchrotron model, and so can not provide the peak 
in the optical/NIR that is required to explain many of the 

observed SEDs. 

A number of authors ( ganders et al. 1989) , for example) 



have argued for the existence of a near-IR bump, somewhere 
around 3/im, corresp ondin g to blackbody emission from hot 
dust. Sanders et al. ( |1989D , also mention the presence of a 
local minimum at 1/im, which they note for its "universal- 
ity" . They interpret this to be due to the finite sublimation 
temperature of dust, causing a drop in the blackbody emis- 
sion, combined with the rise in flux of thermal emission from 
hot (T > 10, 000K) gas. 

Any dip observed in our data occurs at wavelengths 
which are too short to be attributed to hot dust. To observe a 
1/im dip, we would need photometry at longer wavelengths. 



9 EFFECT OF EMISSION LINES ON 
PHOTOMETRY. 

As mentioned in Section ^, the presence of a strong emission 
line in the wavelength range of one of the filters will raise the 
SED level above that of the continuum, which will increase 
the systematic error in the photometry. 

To test this effect, and to see how the fits to the photom- 



etry are affected, we examined several sources whose spectra 
showed significant emission lines. To calculate the contribu- 
tion of the emission line, we define the total flux in the line 
as F\i ne and the total continuum flux under the line as F con t. 
Then the change in magnitude is given by 



Am = 2.5 log^ 



2.5 log ! 



Fcont 
AA + Wy, 

A\ 



where W\ is the equivalent width of the line, and AA is 
the wavelength range over which the flux is measured. The 
equivalent widths and fluxes were computed using the splot 
routine in iraf. 

Note that we are only able to do this analysis for a small 
number of sources, as we do not have spectra for all sources, 
and those spectra that we do have are often poor quality, 
and in almost all cases non-simultaneous. 

The values of the change in magnitude, when considered 
over the same wavelength range as that of the broad-band 
filters (AA ~ 1000 - 2400A for the optical bands), ranged 
from 0.05 - 0.25 magnitudes. By artificially reducing the flux 
in the relevant band by this amount, we could evaluate the 
change in % 2 caused by the presence of the emission line. 
This was done for several sources that had strong emission 
lines present in their spectra (the spectra for a large number 
of PHFS sources will be presented elsewhere (Francis et al., 
in preparation)). A few specific examples: 

• 1510—089 has a strong Ha line in the I band (Am = 
0.14), as well as a combination of a strong H/3 line and promi- 
nent Fen emission in the R band (Am = 0.11). The removal 
of flux corresponding to these lines caused a reduction in 
X 2 jv of nearly 50% (1.28 - 0.68), without greatly changing 
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the resultant fit, although the fitted synchrotron peak was 
at a slightly longer wavelength. 

• 1725+044, similarly, has Ha in / band (Am = 0.13) and 
H/3 + [Oin] in R band (Am = 0.07), and removal of these 
more than halves x 2 l v from 2.49 to 1.17, without changing 
the location of the peak of the synchrotron. 

• 1036—154 has a large Mgn line in the B band (Am = 
0.17), which causes a noticeable upturn in the SED. Re- 
moval of this line reduces the X* l v from 0.77 to 0.56, with- 
out changing the location of the peak wavelength, although 
the 0.5/xm ratio increases slightly (from 5.3 to 7.2). 

• 1136—135 is initially fit with a pure power law, and this 
remains the case after removal of the Mgn line (Am = 0.09 
in B) and both the H/3 and the strong [Oin] (total Am = 
0.11 in /). The power law index softens slightly (from -1.85 
to -1.81) and the yf/v value decreases from 0.43 to 0.29. 

In conclusion, by selecting quasars with strong emission 
lines, we have demonstrated that the largest changes to the 
SED are Am ~ 0.25 in one wave-band. In most cases investi- 
gated, removal of the line flux improved the X 2 l v value, but 
did not significantly alter the nature of the fit. When good 
quality long-wavelength spectra (preferably at least quasi- 
simultaneous with the photometry, which is not the case 
here) are available for these quasars, it should be possible 
to recalculate the fits, taking the emission line contributions 
into account, although we do not expect the general conclu- 
sions to change. 



10 DISCUSSION 
10.1 Types of sources 

The fitting of a single power law to the photometry, as de- 
tailed in Section ^ separates the sources into classes depend- 
ing on the goodness of the power law fit and the value of the 
spectral index. The combined model fits have refined this 
description so that we are able to talk about the different 
types of sources present in the sample. 

Firstly, the bluest sources are all fitted with the power 
law model, and are generally consistent with being the same 
type of sources as optically-selected quasars. The power law 
model also best fits a number of objects of intermediate 
slope (a ~ —1.2), as well as a few red so urce s (a > —0.5). 
These latter sources, 



7.1 



are generally 



as seen in Section 
high polarisation sources, and so are likely to be synchrotron 
dominated (they of course show no evidence for a turn-over 
similar to that being fit by the synchrotron model, although 
their slope could be due to a very steep energy distribution 
(i.e. large p) or a more gradual turn-over). 

The presence of the intermediate power law slopes raises 
interesting issues. Do we see similar sources in optical quasar 
surveys? In Fig. [ll], we have plotted the power law index dis- 
tributions for the sources best fit by the pow er law model, 
as well as the LBQS quasars from Francis (1996). While 



we cannot say that the two data sets come from different 
parent distributions (the Kolmogorov-Smirnov test proba- 
bility for the null hypothesis that the parent distributions 
are the same is 11.3%), there does appear to be an excess of 
quasars in our sample for —1.4 a ^ —0.8. This may indi- 
cate that the accretion disk emission in radio loud quasars 



Power low fits 
LBQS Quasars 





Figure 11. Histograms of power law indices for the power law 
source s in this paper and for the LBQS quasars from Francis 
(Il996|). 



has a broader range of colours, which would have implica- 
tions for models of accretion disk emission. 

The power law sources, like all sources in the PHFS, 
have a radio synchrotron component. The peak of this com- 
ponent may be at longer wavelengths than those covered by 
our observations, i.e. in the IR. This would be consistent 
with obse rvations of flat-s pectrum radio quasars in the far- 
infrared ( Haas et al. 1998 ). These longer wavelength peaks 
would be due to lower energy emitting particles, and would 
therefore constitute the extension of the energy distribution, 
of which we see the high energy end in the two-component 
sources. To test this assertion, one needs to obtain observa- 
tions at longer wavelengths than the K band. We are in the 
process of obtaining observations of a number of sources at 
L band, and the results will be published and discussed at 
a later date. 

Secondly, the combined model generally best fits the 
redder sources, although some intermediate sources are also 
fit, particularly those with SEDs that show inflections. The 
synchrotron components fit to these sources' photometry all 
peak in the rest-frame optical and NIR. 

Thirdly, the other category of sources are those that are 
optically faint. These sources all decrease in flux towards 
shorter wavelengths (indeed, some are not detected in the B 
band), and we suggest that these sources are heavily dust- 
reddened. We note that an exponentially reddened power 
law has the same form as the synchrotron component we 
are fitting, and these sources are fit by a combined model 
with a dominant synchrotron component. 

Finally, we can check the consistency of our model fits 
by plotting the sources on a colour-colour diagram. In Fig. |l^ 
we plot the J — K colours of all sources against the B — 
I colours - that is, the infrared colour against the optical 
colour. We separate the sources into their fitted model types, 
separating the power law sources by their fitted slope (using 
q = —1.3 as the dividing line). 

Clear distinctions can be made between the different 
model types. The blue power law sources lie at the bottom 
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Figure 12. Infrared (J — K) colour versus optical (B — I) colour for all quasars with good model fits. The best fit model is indicated 
by the symbol used. 



left corner, indicating blue colours in both optical and in- 
frared, while as the power law becomes redder, the sources 
move towards the upper-right. Many of the combined fits 
are in the optically-red region of the plot, indicating that 
the SED is turning over in the optical (similar to 1256—229 
in Fig. |^, for example). The reddest sources in B — / are 
sources with optical continua that drop towards the blue, 
in the manner of 1706+006 (see Fig. H), which are typical 
of dusty sources, or sources dominated by host galaxy emis- 
sion, and are fit by a dominant synchrotron component. Note 
that 1706+006 is the faintest of the sources shown in Fig. || 
and so is a good candidate for dust extinction. Other com- 
bined sources, however, are among the bluer sources in the 
optical, but have redder NIR colours than the power law 
sources. These are the sources that show an inflection, simi- 
lar to 1546+067 in Fig. |§|, where the synchrotron component 
is dominant in the NIR but turns over and has less effect in 
the optical. 



law SED. We fitted our models to the observed wavelengths 
of these objects, to investigate the bias created by excluding 
them. Note that this fitting does not take into account the 
change in shape of the synchrotron spectrum due to the red- 
shift of the source. The results are summarised in Table EL 
As can be seen, all but three of the sources are fit with 
the combined model, and the three that aren't are fit with 
quite red power laws. The location of the peak wavelengths 
are generally into the NIR, which are longer than the bulk 
of the distribution of the sources with redshifts (see Fig. ^). 
Also, the synchrotron component that is fitted is generally 
fairly dominant, as evidenced by the fraction values F0.5. 
However, it is unlikely that many of these sources would 
truly be synchrotron sources as they show more the char- 
acteristic shape of dust absorption and are quite faint in 
the optical (again, with the exception of 0048—097, which 
was fit with the power law tail of a dominant synchrotron 
model). We propose that these sources are dust-dominated 
rather than synchrotron-dominated. 



10.2 Sources excluded from the sample 

As explained in Section ^, there are a number of sources that 
were excluded from the sample. Some of these were excluded 
for the sole reason that they did not have a measured red- 
shift. These sources were generally quite red and optically 
faint (thus explaining their undetermined redshift), with the 
exception of the BL Lac object 0048—097, which has a power 



10.3 Sources without good fits 

Up to this point, we have only discussed the results for the 
sources with good (i.e. acceptable at the 99% confidence 
level) fits. A total of 24 sources (or 21 per cent of the to- 
tal) are not fit well by either the power law model or the 
combined model. What sort of sources are these? 
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Table 1. Fitting results for sources that were excluded from the sample purely due to their undetermined rcdshift. 
The best fit model is given along with either the peak synchrotron wavelength (for combined fits) or the power law 
index (for power law fits). F0.5 denotes the fraction of the total flux made up by the synchrotron component, at an 
observed wavelength of 0.5/xm. Note that a * denotes a V — K measurement, not B — K, since a B measurement 
was not made (due to the faintness of the source in B). 



Source name 


B (mag) 




B - K 


Best fit model 


oc 


A p (/jm) 


F0.5 


good/bad 


0048-097 


16.12 




3.61 


Combined 




0.15 


0.97 


good 


1110-217 


24.41 




7.67 


Combined 




1.85 


0.88 


good 


1156-094 


21.95 




5.08 


Combined 




1.04 


0.69 


good 


1648+015 


21.87 




5.31 


Combined 




1.17 


0.82 


bad 


1732+094 


B > 23.5, V = 


21.15 


> 7.13 


Power law 


— 0.39 




0.00 


bad 


2056-369 


V > 23.5, R = 


23.45 


> 5.39* 


Power law 


0.59 




0.00 


good 


2245+029 


21.71 




6.22 


Combined 




1.85 


0.41 


good 


2337-334 


22.93 




6.54 


Power law 


0.36 




0.00 


good 


2344-192 


23.52 




6.27 


Combined 




1.65 


0.83 


bad 



A few sources (~ 8) have roughly power law SEDs, but 
with a little curvature ('n'-shaped) in the blue end of the 
optical. This may be indicative of a small amount of dust 
attenuation or extinction. 

Most of the other sources have one or more photometric 
points that do not smoothly connect with the rest of the 
SED. It is possible that, for these point s, a t least one of 
the systematic errors discussed in Section 2.2 is dominating, 
over and above the level we assigned. In some sources, here 
may also be further emission processes present that we have 
not modelled. 



constant) has a much less drastic effect than changing the 
value of p. As the slope becomes steeper in f\ (that is, bluer), 
we see that only a couple of sources fit by the power law 
change to be fit by the combined model. The distribution 
of A p values is not changed considerably by the variation of 
the as values. 

In summary, reasonable changes to the fiducial values 
of as and p make little difference to the fits, both in the X p 
distribution and the numbers fitted by each model, although 
the difference will be accentuated if the values of as and p 
are both taken to extremes. 



10.4 Alternative values of p and olb 



10.5 Nature of synchrotron model 



Throughout this paper, we have used values for the electron 
energy index of p = 2.5, and for the BBB spectral index of 
as = —1.7. We consider here the effect that changing these 
values has on the results. In Table |^, we list the numbers of 
sources best fit by each of the two models, for each set of 
parameters, as well as the total number of sources fit by one 
of the models. 

Reducing the value of p means that the synchrotron 
spectrum has a bluer slope, which enhances the effect of the 
turn-over. Using a lower value of p in the combined model 
results in some sources, otherwise fit by the power law, being 
instead fit by the combined model. This is indeed seen in 
Table ^, where the number of sources fit by a power law 
decreases as you move from p — 3 to p = 2.5 to p = 2 (if the 
value of as is kept constant), while the number of sources 
fit by the combined model increases. Note also that more 
sources have good fits for the lower values of p. 

The A p distributions are also affected by a changing 
p value. For lower p values, there are more sources with 
shorter X p values (that is, close to 0.1/im). This is due to the 
slope of the power law tail of the synchrotron spectrum: for 
higher p values the slope is redder, and so such a synchrotron 
spectrum peaking around O.l^tm would contribute too much 
at the longer wavelengths. 

For the spectral index of the blue optical power law 
used in the combined model, we considered both ste eper 



and flatter values than the one mentioned in Section 4.1 



The dispersion found by Francis ( 1996 ) for the slopes of 
LBQS quasars is ± ~ 0.3, and so we consider here slopes in 
fx of as = —1.4 and as = —2.0. 

Changing the slope of this power law (while keeping p 



Thus far, we have been considering a synchrotron model 
that peaks at some A p and then turns over sharply (i.e. ex- 
ponentially). The reasons we chose this were for consistency 
with other modelling done for optical synchrotron emission 
( Meisenheimer et al. 1996), and to provide the maximum 
contrast with the power law model. However, synchrotron 
emission could alternatively be present in the form of a 
power law. What if we have a model of the same form as 
the combined model, but with this power law synchrotron 
model instead? Is this any better at fitting the observations? 

We constructed such a model, being a linear combina- 
tion of the blue (BBB) power law from earlier, and a power 
law of variable index: /(A) = aA~ 1,7 + b\ as . To distinguish 
it from the BBB power law, we restricted the indices to the 
range —1.6 < as < 1.5. This model was then fitted to the 
photometry in the same way as previously. Note that the 
special case of a = is simply the power law fit from Sec- 
tion |§] (with the index restricted to lie in the above range) . 

The X 2 / l/ values for this model are plotted in Fig. 
against those for the combined model detailed in SectionW 
The symbols used indicate which model is preferred, and 
whether it is a good or bad fit. For sources that are bluer 
than A -1,7 , the best fit is always with no synchrotron compo- 
nent (for both forms of the synchrotron model), as the pres- 
ence of the second component will redden the blue power 
law. Hence, for these sources, both the original combined 
model and the two power law model give the same fit, con- 
sisting of just the A -1,7 (i.e. BBB) power law. These sources 
are indicated by the "Equal fit" symbols on Fig. |l3[ 

The power law synchrotron model obviously does not do 
as well at explaining the optical/NIR SEDs of the quasars 
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Table 2. Number of sources fitted best by each model. Only those fits accepted at the 99% confidence level are 



shown. Values 


of the 


paramet 


?rs used 


are p = 


2.0, p = 


2.5 and p 


= 3.0, 


as well 


as a.Q 


olb = -2.0. 




















V 


2.0 


2.0 


2.0 


2.5 


2.5 


2.5 


3.0 


3.0 


3 




-1.4 


-1.7 


-2.0 


-1.4 


-1.7 


-2.0 


-1.4 


-1.7 


-2.0 


Power law 


46 


44 


45 


49 


48 


48 


53 


52 


52 


Combined 


49 


51 


50 


44 


45 


44 


39 


40 


40 


Total 


95 


95 


95 


93 


93 


92 


92 


92 


92 
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Figure 13. Values of x 2 / u f° r the combined models featuring 
either power law synchrotron models or synchrotron models with 
a turn-over. The symbols relate to which model is preferred, and 
whether the fit is accepted at the 99% confidence limit or not (i.e. 
good or poor). The line shown is the line of equal x 2 l v - 



as the model with a turn-over. Only 22 sources have a 
good power law synchrotron fit that is better than their 
fit from the synchrotron model with a turn-over, and only 
two of these (1034-293 and 2329-415, with indices of 
as ~ —0.21 and as = respectively) are significantly bet- 
ter. These power law synchrotron sources may have turn- 
overs at shorter wavelengths, beyond our currently available 
data. The only way to tell would be to obtain UV photom- 
etry (ideally contemporaneous with the optical/NIR). 

We note that many of the sources from Fig. |E] that arc 
in the cross-hatched histogram (that is, had better single 
power law fits than combined fits) and had a > —1.6, also 
had fits with the two power law model that were better 
than those with the turnover synchrotron model. (This is 
not always the case, however, since going from one power 
law to two decreases the degrees of freedom, so the x 2 / u 
value will be greater.) In most of these cases, the blue power 
law (the A -1,7 component) is either absent or has a very low 
normalisation. 

For most sources, however, the synchrotron model with 
a turn-over gives a much better fit to the observed opti- 
cal/NIR SEDs (many of these sources show large ratios of 
synchrotron to blue power law flux). This result provides 



clear evidence that there is a break or a cutoff present in 
the energy distribution of the radiating particles. 



10.6 Emission lines and BL Lac objects 

We can use the results of the emission line study in See- 
to draw some conclusions about BL Lac objects, 



7.2 



tion 

which are noted for having very small or non-existent emis- 
sion lines. This is usually explained by the presence of a 
bea med synchrotron compo nent from a jet directed towards 
us ( Blandford fc Rees 1978 , for example). The results shown 
in Fig. |1C| seem to support this explanation, as the equiva- 
lent widths tend to decrease with increasing amount of fitted 
synchrotron, as you would expect if there was an extra syn- 
chrotron component over the top of the emission lines. 

Also, the fits to the BL Lacs in our sample are sug- 
gestive of strong synchrotron components. All nine BL Lacs 
are well fit by the combined model, and the 0.5/mi ratios 
of the synchrotron to power law components range from 3 
(0537-441) to over 700 (1256-229) (see Fig. @). These ra- 
tios, on average, are higher than the rest of the population 
(although it is only a small group). 

However, there may be more to the story. As can be 
seen from Fig. |^, there are a number of sources with com- 
parable ratios to the BL Lacs. Also, Fig. [l(] shows signs of a 
large dispersion in the equivalent width distribution, partic- 
ularly at the high-ratio end, with sources present that have 
relatively high ratios at the emission lines in question, as 
well as large equivalent widths. Perhaps there are other fac- 
tors that determine whether an object is a BL Lac or not. 
One factor that is likely to be important is the strength or 
dominance of the emission line region, with BL Lacs having 
an intrinsically weak emission line flux. Such an idea agrees 
with other studies: for example, Ghisellini et al. ( 1993| ), who 
found that BL Lac objects had significantly lower Doppler 
factors than core-dominated quasars. This was interpreted 
as suggesting that BL Lac objects have intrinsically weaker 
emission lines than core-dominated quasars. 



11 CONCLUSIONS 

The results of this paper can be summarised in the following 
conclusions: 

(i) Radio-loud quasars require one of at least two com- 
ponents to fit the rest frame optical SEDs. These two com- 
ponents are well modelled by an optical power law which 
in some sources is similar to the blue power law seen in 
optically-selected quasars, and a synchrotron component 
that turns over in the range 0.1/^m JS \ p & 3fim. 
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(ii) The model fits require the synchrotron component to 
have a break or turn-over at some wavelength A p in the ob- 
served wavelength range. A synchrotron spectrum with this 
break gives a much better fit to the data than a pure power 
law synchrotron model. The value of A p provides strong con- 
straints on the modelling of particle acceleration and emis- 
sion mechanisms. 

(iii) There is more than 4 orders of magnitude variation 
in the ratio of the synchrotron component to the power law 
component in those sources fit by the combined model. All 
the sources identified as BL Lacs lie at the high end of the 
distribution. 

(iv) The fitting procedure is robust to reasonable changes 
in the values of the model parameters. The fits are statisti- 
cally consistent with the data. 

(v) There may be an excess of red power law sources com- 
pared to optically-selected quasars. This will have implica- 
tions for emission models of accretion disks. 

(vi) Sources which are well-fitted by power laws have a 
synchrotron component which may peak in the IR. In the 
optical, these sources may be similar to optically-selected 
quasars. 

(vii) Optical polarisation measurements (taken from the 
literature) support the synchrotron model fits to the SEDs. 

(viii) There is some evidence for the equivalent widths 
of emission lines being reduced by the presence of the syn- 
chrotron component. This supports the hypothesis that BL 
Lac objects have a dominant synchrotron component, al- 
though other effects may contribute to the lack of emission 
lines in BL Lacs. 

(ix) The red component fitted to the data, which has 
a turnover at some peak wavelength, cannot be hot dust. 
Emission from hot dust cannot peak at A rcst < 1.66pim as is 
required by our data. 

(x) The red sources excluded from our analysis may pro- 
vide interesting examples of extreme sources in the sample, 
e.g. very dusty sources, high redshift sources etc. 
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